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Abstract 



In this work we deal with the problem of support estimation under shape restrictions. The shape 
restriction we deal with is an extension of the notion of convexity named a-convexity. Instead of 
assuming, as in the convex case, the existence of a separating hyperplane for each exterior point 
we assume the existence of a separating open ball with radius a. Given an a-convex set S, the a- 
convex hull of independent random points in S is the natural estimator of the set. If a is unknown 
the r„-convex hull of the sample can be considered. We analyze the asymptotic properties of the 
r„ -convex hull estimator in the bidimensional case and obtain the convergence rate for the expected 
distance in measure between the set and the estimator. The geometrical complexity of the estimator 
and its dependence on r„ is also obtained via the analysis of the expected number of vertices of the 
r„ -convex hull. 
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1 Introduction 

Let 5 be a convex set in the plane. Starting from the classical papers by |[T4l [T5l . asymptotical 
behavior of the convex hull of random points in S has received great attention. Also, expressions 
for the expected area, perimeter, and number of vertices of the convex hull of a sample have been 
object of research. From the point of view of set estimation, the convexity assumption has been 
extensively considered in the literature. If we assume that the set of interest S (for instance the 
unknown support of an absolutely continuous distribution) is convex, then the convex hull of a 
sample from that distribution turns out to be a good choice to recover the shape of the support. 
fS J carry out the asymptotic analysis of the convex hull estimator for general dimension (in terms 
of the Hausdoff distance between the estimator and the set). Computations of measures of the 
convex hull, such as the number of vertices or the volume become quite complicated. In fact, most 
of the known results concern the asymptotic behavior of the expected value of some interesting 
geometrical characteristics such as the area, perimeter or the number of vertices. Only recently the 
asymptotic analysis of the variance or the limit law of these quantities have been performed for 
general convex sets and dimension d see, for instance, |[T2l . We refer to the surveys by IITtI for the 
classical results on convex set estimation and fT3l for more recent results on the subject. 
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Convexity can be a restrictive assumption. Just as an example, it would limit the support esti- 
mation problem to connected supports, which is clearly inadequate, for instance if several groups 
are presented in S and we are interested in performing a cluster analysis. On the other hand, using 
the convex hull as an approximation of a non-convex set leads to considerable errors in the estima- 
tion. A milder shape-restriction which appears in set estimation is a-convexity, see 1 19 1. This shape 
restriction assumes that a ball with radius a can roll freely in the complement of S (see next section 
for a formal definition of a-convexity). This work deals with the study of a natural estimator when 
this restriction is imposed, the a-convex hull of the sample, that is, the smallest a-convex set which 
contains the sample. If a is unknown, we may replace a by a sequence of parameters r„ which 
goes to zero as n tends to infinity. Some results about the asymptotic behavior of the r„ -convex 
hull of the sample can be found in |16]. Here, we are concerned with the convergence rate for 
the expected distance in measure between the set and the estimator. We prove that the obtained 
convergence rate is sharp and cannot be improved in general. We also study the dependence on r„ 
of the expected number of vertices of the r„ -convex hull estimator. This quantity provides some 
information about the complexity of the estimator in the sense that the more vertices the estimator 
has, the more complex the estimator is. 

The rest of the paper is organized as follows. The shape restriction and the estimator are 
defined in detail in Section|2] The main results are established in Section[3l All proofs are deferred 
to Section m 

2 The estimator, the shape restriction and general tools 

2. 1 The estimator 

In what follows we assume that 5 is a (nonempty) compact set in the bidimensional Euclidean 
space M?, equipped with the ordinary scalar product (•,•) and norm || • ||. We also assume that a 
random sample of points Xi , . . . ,X„ from a distribution Px with support S is observed. The goal is 
to reconstruct the set of interest S. Several alternatives have been considered in the literature. For 
instance, under no shape restriction on S, [5] and [7] proposed as estimator of S the union of balls 
of radius s„ with centers in the sample points. See IS for some new results about this estimator. 
However, if it assumed that S fulfils some smoothness restriction then a more efficient estimator 
can be provided. Thus, under the assumption that S is convex, the convex hull of the sample is the 
natural estimator. As it was mentioned in the Introduction, this paper focuses on the problem of 
estimating a set under a more flexible assumption than convexity, named a-convexity. A set A is 
said to be a-convex if any point that does not belong to the set is contained in an open ball (not 
necessarily centered in the point) which does not intersect the set. This recalls us the definition of 
convexity and the existence of a separating plane for each exterior point. In fact, a convex set is 
also a-convex for any value of a. From its definition it can be easily seen that a set A is a-convex 
ifA = Ca{A) where 

Ca{A)= n {B{x,oc)y 

{B{x,a): B{x,a)f]A=ld} 
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is the a-convex hull of the set A, that is, the smallest convex set which contains the set. Here B{x, r) 
denotes the open ball with center x and radius r and A' the complement of A. In what follows B 
and B stand for S(0, 1) and 6(0, 1), respectively. Moreover, from now on, A and dA will denote the 
closure and boundary of A, respectively. 

The a-convex hull of a set A can be also written as the closing of the set, that is, 

Ca{A) = {A(BrB)QrB, 

where © and denote the Minkowski addition and subtraction, respectively. For two sets A,C the 
Minkowski addition is defined by A(BC = {a + c : a £ A, c GC} whereas the Minkowski subtraction 
is A e C = {x : {x} © C C A}. For A G M, AC = {Ac : c G C}. See HH for more details on these 
morphological operators. 

Now, let us assume that S is a-convex for some a > 0. Given a random sample =^ = 
{Xi,. . . ,X„} from Px with support S, the a-convex hull of the sample 

Ca{^n) = {^neOtB)QaB 

turns out to be a natural estimator for the set S. This estimator has the drawback of depending 
on the (possibly) unknown parameter a. This difficulty can be overcome by taking a sequence of 
positive numbers {r,,} converging to zero as n tends to infinity. This ensures that rn< cc for n large 
enough and therefore S is also r„-convex. For the sake of simplicity we assume that r„ < a for all 
n and define the estimator 

Sn = CrS^n) = {^n © r„B) © r„B. (1) 

Our goal is to analyze the asymptotic properties of this set estimator. Here we will consider the 
distance in measure to quantify the similarity in content of S and 5,,. As measure we will use the 
Lebesgue measure \l. Hence, the distance between S and Sn is defined as 

d^{S,Sn)=\l{SASn)=\l[(S\S„)\MSn\S))=\l{S\Sn), 

since with probability one ^ C 5, which implies 5„ C S. 

2.2 The shape restriction 

The estimator ([1} was proposed in |[T6l . In that paper the convergence rate for the Hausdorff 
distance is provided, under the assumption that 5 is a smooth a-convex set. Apart from the a- 
convexity of 5, it is also assumed that S'^ is a-convex. Both conditions imply that S belongs to 
Serra's regular model. See |[T9l for an exact geometric characterization of Serra's regular model in 
terms of a-convexity and free rolling conditions. Essentially, a nonempty compact set S belongs to 
Serra's regular model if, for some a > 0, 

(R) A ball of radius a > rolls freely in S and in S'^ . 

We say that a ball aB rolls freely in a closed set A if for each boundary point a G 5A there ex- 
ists some X G A such that a G fi(x, a) C A. Note that the free rolling condition presented here 



3 



is not exactly the same as the one given in [19|. In that paper it is also required that AQ aB is 
path-connected in order to preserve the physical meaning of rolling freely. We have suppressed 
this additional requirement in our definition of free rolling since it will not be necessary for our 
purposes. Condition (R) is enough in order to guarantee that both S and S'^ are a-convex. It also 
guarantees the existence at each point G 55 of a unique outward pointing unit normal vector 77 {s) 
such that 

B{s-ar\{s),a) cS dindB{s + ar\{s),a) cS'-', 

The proof of these geometrical facts, see Appendix A in lITTl . can be thought as an alternative proof 
for Remark 3 in [19] referring to the validity of its Theorem 1 when the set S is not assumed to 
be path-connected. Another implication of Assumption (R) has to do with the concept of positive 
reach of a set, not mentioned in |[T9l . ifTOl defines the reach of a nonempty closed set A in the d 
dimensional Euclidean space, reach (A), as the largest a, possibly infinity, such that if x ^M!^ and 
d{x,A) = inf{ ||x — y\\:y^A}<a, then the metric projection of x onto S is unique. fTOl provides 
a generalization of the Steiner's formula for sets with positive reach. Recall that, roughly speak- 
ing, the Steiner's formula establishes that the Lebesgue measure of the closed r-neighbourhood, 
B(A,r) = {x : d{x,A) < r}, of a convex set A can be expressed as a polynomial of degree at most 
d in r. Federer's result says that the same holds for sets of positive reach and r < reach (A). It can 
be proved that, under Assumption (R), the reach of both 5 and S'^ is greater than or equal to a. 

2.3 Tools: Unavoidable families of sets 

The procedure of bounding the expected value of (i^(5,5„) becomes easier if we replace the 
proposed estimator by 

S„ = {,5^n®rnB)QrnB. (2) 

It is important to note that, although we use the same notation 5„ for both ( ^„ © r„B) r„B and 
( © rnB) Q rnB, both estimators are not necessarily equal, see Figure [T] However, it is not 




Figure \: For the point set ^ = {Xi,X2,X^}, {S: ®rB)QrB = 3t: and {3^^ ®rB)QrB = 3t: VJ{c]. 
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difficult to prove tliat tiiis event lias probability zero, see Appendix B in ifTTl . Hence we can 
compute E((i^(5,5„)) by using either O or Q. Then, we can write 

E{d^iS,Sn)) = Ei^iS\Sn)) = J^Pix^S„)^idx) 

= I^P{3yGB{x,rn):B{y,rn)n^n = ®)^{dx). (3) 

So, the goal is to find a bound for P{3y G B{x,rn) : B{y,r„) n = 0)- This bound will also 
allow us to obtain a bound for the expected number of extreme points of Sn- As in the convex case, 
it is said that a sample point X, is an extreme point if Z, G dSn- The number of extreme points 
provide us with information about the complexity of the estimator. In the convex case, removing 
the extreme points have been used in data depth for ordering multivariate data sets, see [IJ. A 
similar idea can be used in the non-convex case. If J/n denotes the number of extreme points, then 

lE(^) = nP{X„ is an extreme point). 

It can be easily seen that X„ is an extreme point of Sn if and only if X„ belongs to the boundary of 
an open ball with radius r„ which does not intersect So, conditioning on we get 

]E(^) = nf P{3y G 5B(x,r„) : B{y,rn) D ^„_i = (I))Px{dx) 
Js 

< n jp{3yeB{x,rn):B{y,rn)r\^n-i=^)Px{dx), (4) 

where = {Xi, . . . ,X„_i}. Hence, if we were able to obtain an upper bound for P(3j G 

B{x,rn) : B{y,rn) r\3^n=^) we would get a bound both for E(cf^(5,5„)) and E(^). The idea 
for bounding this probability is to make use of the concept of unavoidable family of sets, defined 
below. 

Definition 1. Let x G K^, r > and (fv,r = {^{y^^) '■ y S B{x,r)}. The family of sets ^x,r is said to 
be unavoidable for Sx.r if for all B{y, r) G S'x,r, there exists U G '^x.r such that U C B{y, r). 

As a consequence of Definition [T] if '^x,r„ is a finite unavoidable family of sets for («^,r„, then 
P(3jGB(x,r„):B(3;,r„)n^„ = 0) < P(3[/ G ^,,,„ : n jr„ = 0) (5) 

If we define for each x G 5 a family ^x^r,, unavoidable and finite for then, from (|3]) and (|5]l, 
it follows that 

Eid^iS,Sn))< f £ {\-Px{U)Ypi{dx)< f £ tx^{-nPx{U))pL{dx), (6) 
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where in the last inequahty we have used (l—x) < exp(— x), for < a: < 1. From Q it is apparent 
that the problem of finding an upper bound for E((i^(5,5„)) (the same holds for E(o4f,)) reduces 
to the problem of finding a lower bound for Px{U), for all U G ^x,r„- In view of (O it would be 
desirable that, both the lower bound and the number of elements of the family '^x,r„ > depend in the 
simplest possible way on the point x. In order to find a lower bound for fx (f^ ) it is useful to assume 
that the probability distribution Px is uniformly bounded on S, that is, 

35 >0 such that Px(C) > di^{CnS) 

for all Borel set C C M^. Crearly, this includes the uniform distribution on S. 

3 Main results 

The main theorem of the paper provides the convergence rate of the expected value of d^{S,Sn). 
The concept of unavoidable family, introduced in Section |2j plays a major role in the proof. In 
Theorem ID we show that the obtained convergence rate cannot be improved. 

Theorem 1. Let S be a nonempty compact subset of such that a ball of radius a > rolls 
freely in S and in S'^. Let X be a random variable with probability distribution Px and support S. 
We assume that the probability distribution Px satisfies that there exists 5 > such that Px (C) > 
5/1 (C n 5) for all Borel subset C C M^. Let = {Xi , . . . ,X„} be a random sample from X and let 
{r„} be a sequence of positive numbers which does not depend on the sample such that rn< a. If 
the sequence {r,,} satisfies 

nr' 

lim-^ = oo, (7) 

n^oo logn 

then 

Eid^{S,S„)) = o(^r;~'n--^y (8) 

Remark 1. proves that, if S is under the conditions of Theorem\l\and {r,,} is a sequence of 
positive numbers satisfying ((71), then, for the bidimensional case, d^{S,Sn) = 0{r^^{logn/n)^^^), 
almost surely. The convergence rate ofE{dj^{S,Sn)) obtained here is, therefore, faster than the 
obtained almost sure convergence rate ofd^{S,S„). Note that the logarithmic term vanishes in 
Moreover, the penalty factor rn is asymptotically smaller than r,7'. 

Remark 2. The proof of Theorem\l\relies on Propositions\J\and\2\(see Section^ which provide 
suitable unavoidable families of sets both for points far from the boundary ofS and close to it. Most 
of the results can be easily extended to the general d-dimensional case. However, some proofs are 
much more involved and of less geometrical nature. The main difficulty in analyzing the general 
case is in proving Proposition^ see [11 /. 

Next theorem shows that the rate in Theorem HJcannot be improved. 
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Theorem 2. Under the conditions ofTheorem\l\ there exist sets Sfor which 

liminfrl nm{duiS,Sn)) > 0. 

Finally we provide a bound for the expected number of extreme points. Note that the bound in 
^ for the number of vertices is almost the same as the quantity which is bounded for the distance 
in measure, see ([3]). The main difference is that (|4]i involves an integral with respect to Px whereas 
([3]) involves an integral with respect to /x. In order to bound integrals with respect to Px by integrals 
with respect to ^ we assume that Px also satisfies 

3j8>0, suchthatPx(C) < i3At(Cn5). 

Again the uniform distribution satisfies the above assumption. 

Theorem 3. Let us assume that the support S and the sequence r^ are under the conditions of 
Theorem [7] Let assume us that the probability distribution which generates the sample satisfies 
that there exit 5,j8 > such that 5lJ.{SnC) < Px{C) < ^lx{Sr\C). Then, 

E(^)=C»(r;^«5). 

4 Proofs 

Theorem\l\ As it was mentioned in Remark |2l Theorem [T]relies on Propositions [T] and |2l Proposi- 
tion [T] gives the desired unavoidable families for the points which are far away from the boundary 
of S. By points which are far away from the boundary we mean those points jc G 5 such that 
d{x,dS) > r„/2. Taking into account Definition[T] it will not be difficult to define a suitable family 
'^x,r„ in this case. We need that, given y £ B{x,r„), there exists U G '^x,r„ such that U C B{y,r„). 
It would be also desirable that U was totally contained in S and that IJ.{U) was of the maximum 
posible order rf^. This would ensure the best possible rate for Px{U). Note that if x G 5 and 
d{x,dS) > r„/2, then the ball B{x,r„/2) is fully contained in S. So, the idea is to divide B(x,r„/2) 
into a finite number of subsets. Here, we will consider a partition of B{x, r^/l) into circular sectors. 
The choice of circular sectors rests upon two main reasons. First, the measure of a circular sector 
of B(x, r„/2) is of order r,^. Second, if the central angle of the defined sectors is not too large, then 
the resulting family ^x,r„ is unavoidable. 

Before the statement of Proposition [T] we give the precise definition of the circular sectors and 
introduce some basic notation that will be useful later. Thus, let §2 = {" ^ • II "II = 1} denote 
the unit circle in and = (0, 1) G M^. Let % y be the angle between the (nonzero) vectors u 
and V. It is understood that (pj, „ G [0, n] and (p„ v = (pv,» For m G §2 and 6 G [0, 7r/2], we define the 
cone = {x G : {x, u) > \\x\\ cos 6} and the circular sector C^^ = n B(0, r). Note that Cf^ 
is the circular sector with central angle 26 enclosed by the radii vi = rMg{u) and V2 = r^Q^{u), 
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where !Mq : — > denotes the counter-clockwise rotation of angle Q, whose associated matrix 
with respect to the canonical basis is 

/ COS0 — sin0 \ 
Y sin0 COS0 y ' 

In Figure |2] we show an example of C^^. 





lu/ W V2 







Figure 2: Circular sector C„ ^. 



Proposition 1. Under the conditions of Theorem^ for all x & S such that d{x,dS) > r„/2, there 
exists a finite family ^,r„ with m\ =6 elements, unavoidable for Sx,rn cit^d that satisfies 

Px{U)>Lirl U€^,,r„, 

where the constant Li > is independent ofx. 

Remark 3. This proposition can be easily generalized for dimension d. The main difference is that 
m\ is in general unknown since it depends on the number of cones we need to cover the unit ball. 

Proof. First consider the family ^Q^r„ = {C^^^' " ^ where W (I'M} denotes a set of unit 
vectors that divides the unit circle into six circular sectors with central angle 7r/3. Figure |3] shows 
one possible choice of W and the corresponding family ^o,r„ ■ To simplify notation somewhat, we 
abbreviate Cu^^ and df/-,^ to C„ and Cu,r„ , respectively. Note that the definition of W implies that 

B(0,r„)= U Cu,r.r 

The fact that '^o/,, is unavoidable for <§qj-„ easily follows from Lemma [2l stated below. To see 
this, note that for B(y,rn) G <^o.r„, there exists u ^ W such that y G Cu,r„- Now, by Lemma |2l 
Cu,r„ C B(y,rn) and therefore C„ r„/2 C B{y,r„). This completes the proof that ^o,r„ is unavoidable. 
Thus, it remains to prove Lemma|2] First we establish, without proof, Lemma[T]which characterizes 
the points in and simplifies the proof of Lemma |2] 
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Lemma 1. Let x^Q. Then 

We are now ready to state and prove Lemma |2] This lemma reveals that the partition of B(0,r„) 
into circular sectors with central angle 7r/3 is indeed a sensible choice, since it guarantees that 
^o,r„ is unavoidable. 

Lemma 2. For all m G §2 '^nd r > 0, 

Cu,rC n B{y,r). 
yec,,,- 

Proof. Let z € C„,,- We need to show that, for all y G Ci,,r, \\z — y\\ < r. Assume, without loss of 
generality, that z and y are both non zero vectors since the result is trivial otherwise. We have that 

Ik -yf = llzf + Ibf - 2||z|| Ibll cos 

By the triangle inequality for angles and Lemma[T]we have (p, y < (p^ „ + 9„ y < -j. Hence, 

llz-jf <lkf + lbf-|Nlll3'll<max(||z||Mbf)<r2. 

□ □ 

Once we have proved that r„ is unavoidable for (^o,r„ consider, for each x G 5 such that 
d{x,dS) > rn/2, the family %,r„ = {4©%r„ = {{x} eCu,r„/2, " S W}. The family %^r,„ ob- 
tained by translating the family ^o^r„ by the vector x, is unavoidable for (f^ ,.„ > as we state in Lemma 
[3i We skip the proof since it is straightforward. 
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Figure 4: For x G 5 under the conditions stated in Proposition [7] we have that 
{4ffiQr„/2 CB(x,r„/2) C5. 



Lemma 3. Let be an unavoidable family for (fo,r- Then ^ = {x} © ,- = {{x} (BU, U G 
^o,r} '■5' unavoidable for S'x^r- 

To complete the proof of Proposition [T] it remains to give a lower bound for the probability of 
the sets of the unavoidable family we have just defined. For each m G ^ we have that 



This follows simply because {x} ©C„ r„/2 C B(x, r„/2) C 5 since <i(x, dS) > r„/2 and the Lebesgue 
measure is invariant under translations, see Figure |4] Therefore, 



Before proceeding to the definition of unavoidable families of sets for points x G 5 with d{x, dS) < 
r„/2, we wish to emphasize some aspects of this kind of families. Recall that for points which lie 
far away from the boundary we have proved that it is enough to consider circular sectors with radius 
r„/2 and central angle 7i/3. Using the same argument for points x£S such that p = d{x, dS) < r„/2 
we only could infer that B{x,p) C S and hence the lower bound for the probability of these circular 
sectors would be of order p^. However we can find larger unavoidable sets and improve this bound. 
To see this, assume without loss of generality that x = and divide B{0,r) into a finite number of 
sectors C^^ with 6 > 0. Then for fixed u, 



PX (W ® Qm-„/2) > 5^ (W ®Q,r„/2 n5) = 5 IX (W ® Q,,r„/2) = (C„,,„/2) . 




for Li = dn/24 > and the proof of Proposition [T]is complete. 



□ 
□ 



U= f] B{y,r) 



(9) 
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is the largest set contained in B{y, r) for all y G Cf ,.. The measure of U depends on 6. For example, 
if 6 = 7i/2 then we divide B{0,r) into two circular sectors with central angle 7i. In that case, it can 
be easily proved that U = {0}. Smaller values of 6 result in larger sets U. In particular. Lemma 
[2] shows that, fixed 6 = n/6, the set in ^ contains at least one circular sector with central angle 
n/3. In Proposition |2] we show that for points x £ S with p = d{x, dS) < r„/2 and 6 = 7i/6 we can 
give a lower bound for Px{U) of order r//^p^/^. Note that this bound is better than the one we can 
obtain for circular sectors of 5(x,p). Hence, Proposition [2] provides the second key result in the 
proof of Theorem [T] At this point it is worth discussing some of the properties of the sets 

P B{y,r), with u G §2, and r > 0. (10) 

As we show in Lemma |4] below, these sets are known in the literature as Reuleaux triangle, see 
Figure[5] They solve the problem of finding unavoidable families of large sets for the bidimensional 
case. One can be tempted to generalize the idea for the (i-dimensional case. However the argument 
in is somewhat different since it becomes tough to handle with the intersection in ([TOl l when 
d > 2. Note that it is fundamental not only to define large unavoidable sets but also to measure 
them. This causes technical difficulties as the dimension increases. 



u / 




Figure 5: Reuleaux triangle. 



Lemma 4. Given u G §2. we have 

n B{y,r)=B{0,r)nB{vi,r)nB{v2,r), 

yeCu,r 

where vi = r^{u) and V2 = r^^^{u), ^ : — > M.^ being the counter-clockwise rotation of angle 
7r/6. 

Remark 4. As previously discussed, the set B(0,r) nB(vi,r) nB{v2,r) in M.^ is the so-called 
Reuleaux triangle. Formally, the Reuleaux triangle is defined from an equilateral triangle with 
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sides of length I. It is constructed by drawing the arcs from each polygon vertex of the equilateral 
triangle between the other two vertices. Thus, the Reuleaux triangle is the set bounded by these 
three arcs. An important property is that it is a set of constant width I, see Figure It is known 
that the diameter of a set of constant width I is precisely I. See Ml/. /El/, (md the references 
cited therein for a detailed development of these concepts. 

Proof. It is straightforward to verify 

n B{y,r)^B{0,r)r\B{vur)r\B{v2,r). (11) 

Let us now consider the reverse content. Let x G B(0, r) n B{v\ , r) r\B{v2,r) and y G Cu,r- We need 
to show that Hjc — j|| <r. It follows from (fTTI ) that y G B{0,r) nB{vi,r)nB{v2,r) and hence, since 
the diameter of the Reuleaux triangle is r, the result holds. □ 

□ 



'I o ^ 

Figure 6: Sets of constant width. 

We now concentrate on the points x which are close to the boundary of 5. Recall that by points 
which are close to the boundary of S we mean those x ^ S such that d{x, dS) <r„/2. As previously 
described, we shall consider in this context unavoidable sets which are larger than the circular 
sectors used for points away from dS. The unavoidable sets U we shortly define guarantee a lower 
bound for Px{U) of order rn^d{x,dS)^^^. Proposition |2] makes these ideas precise. 

Proposition 2. Under the conditions of Theorem\l\ for all x ^ S such that d(x,dS) < r„/2, there 
exists a finite family ^,r„ with m2 = 6 elements, unavoidable for Sx,rn <^f^'^ t^'^t satisfies 

Px{U)> L2rl d{x, dsf^, Ue '^,,r„ , 
where the constant L2 > is independent ofx. 

Proof. Let x G 5 such that p = d{x,dS) < r„/2 < a. Since reach (y) > a there exists a unique 
point Ppx ^ dS such that p = \\x — Prx\\. The rolling condition ensures the existence of an unique 
unit vector T] = 77 {Ppx) such that B{Prx — arj ,a) C S and therefore, given an unavoidable family 

Px{U)>5ix{UnS)>5n{Ur\B{Prx-ari,a)), U e%,r„- (12) 
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Note that this simpUfies the proof since by ([T2] | it follows that we just need to define a suitable 
family '^^,.„ and bound /x(?7 H B {Prx — ari,a)) for U G'^v,r„- Let us consider a composite function 
T formed by first applying an orthogonal transformation ^ : — > such that ^(^2) = — ^ and 
then applying the translation by the vector x, see FigurelV] In particular T{0) = x, T{{a — p)e2) = 
x — {a- p)r\ = Ppx — at], and 

T{B{{a-p)e2,a)) =B{Prx- ar\,a). 



B{{a-p)e2,a) 





Figure 7: For the function T, T {B{{a — p)e2,a)) = B{Prx — arj ,a). 

It can be easily seen that the following result holds. 

Lemma 5. Let ^o,r t>e an unavoidable family for Sqj and let G : — > he an orthogonal 
transformation. Then {^(?7), IJ £ cil^o an unavoidable family for Sq,- 

What Lemma[5] asserts is that the orthogonal transformation of an unavoidable family for (fo,r„ 
results in another unavoidable family for (fo,r„ ■ On the other hand, Lemma |3] established that the 
result of the translation of an unavoidable family for Jo^r„ by the vector x is an unavoidable family 
for Sx,r„- As an immediate consequence, we obtain that ^x,rn = {T{U),U G ^o,r„} is unavoidable 
for (Cr,r„ ■ Furthermore, 

H{T{U)nB{Prx-ari,a)) = n{UnB{{a- p)e2,a)), 

as the Lebesgue measure is invariant under translations and orthogonal transformations. Thus, the 
problem reduces to defining an unavoidable family '%.,„ for (#o.r„ and finding a lower bound for 
H{UnB{{a-p)e2,a)) for all U £ %,r„- 

Before continuing the proof of Proposition |2] it may be useful to make some comments con- 
cerning the measure of the sets U nB{{a — p)e2,cc). Note that when defining unavoidable sets 
for (fo,r„, the main difficulty in giving a lower bound for ^{U nB((a — p)e2, a)) arises with those 
points which lie far away in the direction of the vector — ^2- In fact, 

min n{B(y,rn)nB{{a- p)e2,a)) = n{B{-rne2,rn)r\B{{a- p)e2,a)) 

yeBio,r„) 
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since y = —r„e2 represents the point where the distance between the centres of both balls attains its 
maximum and, as a direct consequence, the intersection its minimum. Recall that, by the definition 
of unavoidable family, for each y e B{0,r„) there exists U G ^o,r„ such that U C B(y,r„). So, it 
is more involved to find unavoidable sets U with large enough IJ,{U nS((a — p)e2, oc)) for points 
close to — r„e2- This motivates dividing B(0,r„) into two subsets as follows B{0,r„) = U=^r„ 
where 

?^,„ = |yGB(0,r„): {y ^e2) > -^-\\y\\}^ and = G B(0,r„) : (y^e2) < -\\\y\}^. 

Figure [8] shows the sets §^r„ and .!fr„- Roughly speaking, contains the points y G B{0,r„) 
for which B{y,rn) nB{{a — p)e2,cc) is small. Therefore, the unavoidable sets U in this case 
should be carefully selected. On the contrary, %^ contains the points y G B{0,r„) for which 
S{y,r„)r\B{{a — p)e2,C)c) is larger. For these points the sets U can be circular sectors. Propo- 
sition [3] shows that IJ.{U r\B{{a — p)e2, a)) is then large enough. 




Proposition 3. There exists a finite set of unit vectors W'^ C §2 with = 4 elements such that, 
for all y G l^r„. there exists u G such that y G C„,r„ C B{y, r„) and 

At (C„,,-„ n B((a - p )e2, a)) > L^rl p 5 , 

where > is a constant. 

Proof Let us consider the set = {(1,0), (-1,0), (1/2, v'3/2), (-1/2, ^/3/2)}. It is straight- 
forward to verify, see Figure |9] that '^r„ = [JueW^ ^u,r„ ■ Therefore, for all y G there exists 
u G W"' such that y G C„,r„- By Lemma|2]it follows that C„ ,■„ C B{y, r„). It remains to find a lower 
bound for Cu,r„ nB((a — p)e2, cx) for u G W'^ . Note that at least half of the set C„,r„ is contained 
in the halfplane Hq = {x = {xi ,X2) G : X2 > 0} and hence it is sufficient for our purposes to 
concentrate on Cu,r„ n//o. 
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Figure 9: Unit vectors = {(1,0), (-1,0), (1/2, ^3/2), (-1/2, ^3/2)} and Cu,r„, far u e W'^ . 



Let V = Y^p(2a — p). By the Pythagorean theorem, it is straightforward to see that v represents 
the distance to the origin from the points such that dB{{a — p)e2,cc) intersects the axis OX, see 
Figure \M It is also easy to show that B{0, v) HT/q C B{{a - p)e2,a). Therefore, for u G W'^ , 
Cw,T„ n//o C Cu^r„ nB{{a — p)e2,cc) where t„ = min(v,r„). This yields, 

^i{Cu,r„nB{{a-p)e2,a))>^{Cu.r.nHo) > ^At(C„,T„) = ^t,? > ^rlJ^p^^ 

This completes the proof of Proposition |3j with = 7r/12 > constant. □ 

□ 

In view of Proposition [3] we define the family '^q^ = {Cu,r„,ii S W^}, formed by = 4 
elements. We now turn to the points in J^r,,- The aim is to define for those points a finite family 
'^0%' such that, for all y G ,^r„, there exists U G ^o^_ that satisfies U C B(y, r„) and 

^{UnB{{a-p)e2,a))>L^rlp'2, yu G%%^. (13) 

At this point, it may be useful to make some comments concerning the main differences between 
'^r„ and ■ One might be tempted to proceed as before for ^i-„ and define the set of unit vectors 
W-^ = {(-1/2,-^/3/2), (1/2,-^/3/2)}. Again we would have that, see Figure [ID (a), = 

If we repeat the sketch of the proof for and define U to be the circular sectors C„_r„ for 
u e W'^ , we could no longer guarantee the lower bound in (|T3) . Note that the intersection C„ H 
B{{a — p)e2,cx) for u G W'^ is considerably smaller than for u G In fact, it can be easily 
proved that, for u G W'^, iU(C„,r„ r\B{{a - p)e2,a)) < V^p^, as it is shown in Figure [TTI (b). 
Therefore, we need to consider different sets U. We have previously discussed the possibility of 
defining unavoidable sets, larger than circular sectors. For a fixed unit vector u, 

U= f] B{y,rn) (14) 

y^Cu,r„ 
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is the largest set such that U C B{y, rn) for all y E C„_r„- Figure [12] shows C„.r„, for an m G W'^ and 
the corresponding set U defined in (fT4l) . Observe that U r\B{{a — p)e2,cc) is clearly larger than 
Cu,r„ n B((a — p)e25 o;). The difference between both intersections will play a fundamental role in 
obtaining the lower bound in $13]) . In fact, it is not necessary to consider the whole U as defined 
in (fT4] ). For our purposes it is sufficient to measure a portion of U nB((a — p)e2,cc). We shall 
consider sets as the one represented in gray in Figure [13] Its measure is large enough to satisfy 
([T3] ). We give the precise definition of this kind of sets in Proposition [4] This solves the problem 
for the points in . 

Proposition 4. There exists a finite family of sets '^Q^n ^^^^ ~ elements such that, for all 
y G ^r„> there exists U G '^^o^^ such that U C B{y, r„) and 

^{UnB{{a-p)e2,a))>L^rip'2, 

with > a constant. 

Proof. First, let us consider the set B{{a — p)e2, a) r\B{—r„e2,r„), which corresponds to the in- 
tersection between two balls of radii a and r„, respectively, being a + r„ — p the distance between 
their centres, see Figure [14] (a). The values of hi, h2 and A in Figure [14] (b) can be deduced from 
the Pythagorean theorem. They satisfy the following equations 

(r„-/ji)2 + A2 = r2, 
(a-/j2)2 + A2 = a^, 
hi+h2 = p. 
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(a) 



(b) 



B{{a-p)e2,a) 




Cu,r„ nB{{a — p)e2,cc) is contained in the rectangle of height p and base \^p. 



By solving the system, 



p{2a-p) p{2r„-p) f— - 

hi = — — ■ -, h2 = — — ■ ,and X = \ lrnhi-h\. 

2(a + r„-p) 2(a + r„-p) V 

We now define the set 

'^{hi) = {xeM^ : -hi<{x,e2) <0]r\B{-r„e2,rn). (15) 

Lemma|6]provides a lower bound for the measure of ^(/ii). 
Lemma 6. Given the previous set ^{h\), then 

Hinh))>^rlpi. 



(16) 



Proof. We have that 

Ai(^(/ii))= r'2yjlrr,y-y^dy. 
Jo 

For y G [0,hi] we have that y < rn, since by construction h\ < p and by assumption p < r„/2. 
Hence, 2r„j — y^> r„y and 

MWi)) > r 2^dy = \r],h\. 
Jo 3 

Moreover, hi > p /2, since r^ < cc and this completes the proof. □ 

□ 
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Remark 5. Note that the exact value of the integral in di6D can be explicitly computed since it 
coincides with the area of the circular segment defined by the chord that joins the intersection 
points ofB{{a — p)e2, (x)nB{—r„e2,r„). Thus, 

/^('^(/Ji)) =r^arccos ^ ^' ^ - {rn-hy)^lrnhi-h\. 

So, we have defined the set "^(/ji), whose measure verifies the statement of Proposition |4] Next 
lemma shows that "^(/ii) is contained in B{{a — p)e2, a). 

Lemma 7. 

^(/ii) cB((a-p)e2,a). 

Proof LetxG'^(/ji). 

||x- (a -p)e2|P = \\x\\^ + {a- p)^ -2{a - p) {x,e2) . (17) 

By definition, x £ B{—rne2,rn) and therefore < — 2r„(x,e2)- Furthermore, by definition, 
{x,e2) > —hi. Turning to ([T7J we get 

\\x- {a - p)e2\\^ < 2r„hi+{a- pf' + 2{a- p)hi 
= p(2a -p) + (a-p)^ = a^. 

□ 
□ 
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Figure 13: The dashed area corresponds to U = riyec,,, ^(y^^n) with u = (1/2, — \/3/2). 

It follows from Lemmas [6] and |7] that 

At(^(/ji)nB((a-p)e2,a)) >Lripi (18) 

In order to complete the proof, it remains to define the family '^o^, mentioned in the statement of 
Proposition m In view of ([TSl l. it seems natural to divide 'io{h\). We denote Qi = {x = (xi^xj) G 
: ;ci >0} and Q2 = {x = {xi,X2) eM.^ : xi < 0}. Then, = (Gi n ^,-„) U n and, in 
the same manner, ^(/ji) = (2i n'^(/ii))U (g2n'^(/ii)). 

Lemma 8. For all y ^ QiH ^r„ we have that 

a-n^(/ji)cB(y,r„), /=1,2. 

Proof. Let G Qif^'^(hi). First, it can be easily proved that Qi n=^r„ =C„.r„5 with m = (1/2, —1/3/2). 
What we need to prove is ;c G Clyec,,,- ^(y^^n)- It follows from Lemma|4]that 

f] B{y,r„) = B{0,rn)nB{vi,r„)nB{v2,rn), 

where vi = rn^{u) = rn (\/3/2, — 1/2) and V2 = rn^^^{u) = — r„e2. We have by definition that 
X G B(y2,r„). Moreover, \\x\\^ <X'^ + h\ = 2r„h\ < r^, since h\ < p < rn/1. Note that the last 
inequality justifies the choice of p < r„/2. And, 
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since < < A < ^r„/2 and —h\ <X2<0, where h\ < p < r„/2. Thus, we have shown 
that X G B(0,r„) nB(vi,r„) nB(v2,r„) and the lemma is proved for 2i □'^(/zi). The proof for 
Q2 n '^{h\) is analogous. □ 

□ 

In view of the previous results we define the family '^q^ = {2in'^(/zi), / = 1,2}, formed by 
= 2 elements. It follows from Lemma [8] that, for all y E ^r„, there exists / G {1,2} such that 
Qi n '^{hi) C B{y, r„). Moreover, by Lemma|6] 

1 3 2 
Lrl p I < At (^(/ii ) ) = £ At (a- n "^(/Ji ) ) . 

!=1 

The symmetry of the set 'to (hi) with respect to the axis OY implies that the orthogonal transforma- 
tion ff:R^ — >R^ such that iff{x) = G(x\,X2) = i-xuX2) transforms 2i n'^(/zi) into Qiri'^ih) 
and then both sets measure the same, that is, 

At(ein^(/ji)) = At(G2n^(/ji)) = ^MWi))- 

By Lemma|7]we further have that, for / = 1,2, 2, n'^(/ji) C "^(/ii) C S((a — p)e2, a) and hence 

M (2; n '^(/z 1 ) n B( (a - p )e2 , a) ) = M (2,- n "^(/z 1 ) ) > L-^r„' p 5 , where = V2/6. This completes 
the proof of Proposition |4l □ 

□ 
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Now, we define ^o,r„ = '^o^„ ^^r„ ■ mentioned at the beginning of Proposition |2l 

'^x,r„ =\J^\ U G ^o,r„} is a finite family with m2 = + = 6 elements satisfying that, for 
each U G %^r„, 



PxiT{U)) > 5iJ.{T{U)r\B{Prx-aT],a)) = 5}x{U r^B{{a- p)e2,a))>L2rlp 
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where L2 = 5min(L^,L'^). This completes the proof of Proposition |2] □ 

□ 

We are know in position to complete the proof of Theorem [T] Recall that, if we define for each 
S n family '^v,r„ unavoidable and finite for Sx,r„ , then 

Eid^iS,S,))< f £ cxp{-nPx{U))nidx). 

We divide S into two subsets S= {x£S: d{x,dS) > ^] U {x £ S : d{x,dS) < '-f] and then 
K{d^{S,Sn)) < / ^ £ cxp{-nPxiU))^idx) 

J{xeS: d{x,dS)>^}ui^^,^^_ 

+ [ y cxp{-nPx{U))^{dx). (19) 

J{xeS: d{x,dS)<'^}^^^^^^ 

For those x such that d{x,dS) > r„/2 we make use of the families '^x,r„ given in Proposition 
[U which ensures the existence of suitable finite families ^/i^x,r„ and provides a lower bound on the 
probability of the sets U, independent of x. Thus, 

/ y tM-nPx{U)Mdx) 

J{xeS: d(x,dS)>'-f]ui^/^^_^ 

< / my&xp{-nUrl)^{dx) = O (&-^'"'''\ . (20) 

J{xeS: d{x.dS)>'^] ^ ) 

For those a' G 5 such that d{x^ dS) < r„/2, we may consider the unavoidable families ^x.r„ given in 
Proposition |2l We have that 

L y exp{-nPxm^{dx) 

J{xeS: d(x,dS)<^}ui^^,^_^ 
< I m2 exp ( —L2nr}d{x, dS) ~- ) }x{dx) 

J{xeS: d[x,dS)<'-^} V / 

g{^(x))pi{dx), 

■^-'([0,r„/2]) 
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where =^ : 5 — )• M is defined as 3^{x) = d{x,dS) and g{z) = m2exp{—L2nr^z-)- It follows from 
tlie cliange of variables formula (see Theorem 16.12 of [4]) that 

/ g{^{x))n{dx)= [ g{p)^^-\dp) (21) 

>/,^-'([0,r„/2]) J[0.rj2] 

where p = and jU=^"^ is the measure on R defined by ij.£^^\A) = /x(5'"^(A)), for A C M. 
The measure jj.-^^^ is characterized by F{z) = IJ.{x G S : d(x,dS) < z}. Since reach(55') > a, 
F{z) is a polynomial of degree at most 2 in z, < z < a, see 1 10|. Therefore, it is a differentiable 
function and F'{z) is bounded on compact sets. In short, we obtain 

g{p)piSr-\dp) = I m2&x^\-L2nr],pi]F'{p)dp 

,r„l2\ ^[0,r„/2] V / 



< K I m2exp -L2nrnP- }dp 







^2" „2 



^ / ' " fnij-^rn Se 'idv = 0[ rn 3 

2^2 



where we have used the change of variables formula v = L2?ir,| p ^ and also the fact that e ^'v 3 dv < 
oo. Turning to the computation of E,{d^{S,Sn)) in ( fT9l ). it follows from (l20l) and (|22]) that 



E(^/^(5,5„)) = O j^e-^'-" +r„ '^-^J . (22) 

Since r„ is bounded by a and ?ir^/logn goes to infinity, we have e^^'"''" = o{r„ ^^^n^^^^). There- 
fore, E,{d^{S,Sn)) = 0{r„ ^^^n^^/^), which completes the proof of Theorem[T] □ □ 

Theorem^ Let S = B(0, a) and assume that the distribution Px is uniform on S. Our aim is to find 
a lower bound for ¥.{d^{S,Sn))- Thus, 

E(J^(5,5„)) = J^P{3yeB{x,r„):B{y,r„)n:^„ = <l))^{dx) 

> i , ^(33^ G B(x, r„) : ^(3;, r„) n ^„ = 0)^1 (t/.^) . 

J{;ie5: d{x,dS)<'^} 

For each x G 5 such that d{x, dS) < r„/2 let rj = x/ \\x\\ and x = (||x|| + r„)T], see Figure [15] A 
simple geometric argument shows that Px{B{x, < 1/2. Since x G B{x, r„) we have 

E(J^(S,5„)) > / fl-Px(B(x,r„)))'V(Jx) 

j{;feS: rf(x55)<f} 

^ / ( -nPx{B{x,rn))\ 

> exp -^--7 ^ ;U((ix) 

J{xeS:rf{x,3S)<r„/2} V 1 " Px{B{x, rn)) ) 

> I exp{-2nPx{B{x,rn)))n{dx). (23) 

J{xeS:d{x,dS)<r„/2} 
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Figure 15: Given x € B(0, a) such that d{x, dS) < r„/2, we define x. 



Above we used the fact that (1 — z)" > exp(— ?iz/(l — z)) for z € [0, 1). In view of (1231 ) we need 
again an upper bound for Px{B{x,r„)). The previous bound Px{B{x,rn)) < 1/2 is too rough for 
our purposes and obviously it can be sharpened. Let us consider the composed function formed 
by first applying an orthogonal transformation iff :M? ^M? such that ^(t]) = —ei and then ap- 
plying the translation by the vector (a — d{x,dS))e2, see Figure [T6l Since the Lebesgue mea- 
sure is invariant under orthogonal transformations and translations we have that /x(B(x, r„) n^) = 
IJ.{B{-rne2,rn)riB{{a -d{x,dS))e2,a)). The set B{-r„e2,r„)nB{{a-d{x,dS))e2,a) is the in- 
tersection of two balls with radius r„ and a such that the distance between their centres is equal 
to a + r,, — d{x,dS). Recall that this set appeared in Proposition |4] Following the notation used 
previously, B{—rne2,r„)r]B{{a — d{x,dS))e2,(x) = '^{hi)LI £/{h2), where 'W{h\) is given by ([T5]) 
and 

i2/(/j2) = {z G : -{hi+h2) < {z,e2) < -hi} nB{{a - d{x,dS))e2,a). 

Recall that the values of hi and /i2 were easily deduced from the Pythagorean theorem by solving 
the system 

{rn-hif + ?i^ = rl 
{a-h2)^ + X^ = a}, 
/ji +/i2 = d{x,dS). 

Thus, 

^ _ d{x,dS)(2a — d{x,dS)) ^ _ d{x,dS)(2rn — d{x,dS)) 
^~ 2{a + rn-d{x,dS)) ' ^~ 2(a + r„ - J(x,5S)) 

Since 'io{hi) and iz/ (/j2) are disjoint, up to a zero measure set, we have 

pi{B{-rne2,rn)r\B{{a-d{x,dS))e2,a))=pi{'^{hi))+pL{£/{h2)). (24) 
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Figure 16: (a) B{x,r„)riS. (b) Result of applying an orthogonal transformation 0' :M? such 
that i?{r\) = —€2. (c) Translation by the vector (a — d{x,dS))e2. In black s^{h2) and in gray 



First, in order to find an upper bound in (|24] |. we shall see tiiat (/j2)) < /^('^(/Ji)). It can be 
easily proved that IJ.{£/ {hi)) = IJ.{M){h2)), where 

£/o{h2) = {z G : < {z,e2) < hi} nB{-{a - h2)e2,a). 

As in Lemnia[6l we have /i(£^(/22)) = /q^ l\Jlay — y'^dy. Using the change of variable l = h2—y, 
and taking into account that 2ay — y^ = — {a —y)'^ we can write 

Ai(M)(/i2)) = J^"^' 2^a'^-{a-h2 + l)^dl = ^ j'^' Vs(l)dl, 

Similarly we have = 2 /q ' y/r(l)dl, where r{l) = rl- {rn -hi+ 1)^. Note that r(0) = 

^(O) = and h2<h\. It is easy to show that s{l) < r{l) and therefore 

At(j#(/i2)) = 2 f'" ^/7{^)dl < 2 ^/Hj)dl < At(^(/2i)). 
JO Jo 

Now, if we return to Equation (|24] |. we get 

H{B{x,r„)nS) <2n{'<f{hi)). 
An upper bound for ju('^(/ii)) can be easily found since 

At(^(/Ji)) = J^' ^{2rny-y^)dy < J^" y^dy = O (rhj) = O (rldix,dS)'2 
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where in the last equality we have used hi < d{x,dS). As a consequence, 

Px{Bix,r„))<Lrid{x,dS)'2, 

with L > a constant which does not depend on x. Finally, if we apply the latter bound to ( |23] ). 
then we have that 

E{dJS,S„)) > [ expl -2nLrid{x,dS)'A n{dx) 

J{xeS:d{x,dS)<rj2} \ J 

= I , g{^{x)Mdx), 

i 3 

where : 5 — M is defined as ^{x) = d(x,dS) and g{z) = exp{—2nLriiz-)- We use the same 
change of variables formula, see dZB, withF(z) = 7r(a2-(a-z)2). So 

rr„/2 / 1 

E(J^(5,S„)) > expl-2nLrip-2jF'{p)dp 
= J exp ( -2«Lr„-p2 1 2;r(a -p)Jp > na J exp l-lnLrii j dp. 
A straightforward calculation shows that 

E{dJS,Sn))>Crn'n--i f ^ e^'v-Uv 

Jo 

for some constant C > 0. Since nr^ — )• oo, we have 

liminf 4 nm(da(S,Sn)) >C f e^'v^^dv > 0. 

This completes the proof of Theorem |2] □ □ 

Theorem\3\ Note that under the stated assumptions on Px we have that, for any measurable non- 
negative function in 5, (p, we have that 

q){x)Px{dx) < P J (p{x)}x{dx). 

Using this fact in Q, we can follow the same lines as in the final part of the proof of Theorem [Jto 
easily conclude the result. □ 

□ 
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